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We have measured the heat capacity of an optically-trapped, strongly-interacting 
Fermi gas of atoms. A precise input of energy to the gas is followed by single- 
parameter thermometry, which determines the empirical temperature param- 
eter f of the gas cloud. Our measurements reveal a clear transition in the heat 
capacity. The energy and the spatial profile of the gas are computed using a 
theory of the crossover from Fermi to Bose superfluids at finite temperature. 
The theory calibrates T, yields excellent agreement with the data, and predicts 
the onset of superfluidity at the observed transition point. 

Strongly-interacting, degenerate atomic Fermi gases provide a paradigm for strong in- 
teractions in nature ©. In all strongly interacting Fermi systems, the zero-energy scattering 
length is large compared to the interparticle spacing, producing universal behavior Predic- 
tions of universal interactions and effective field theories in nuclear matter (T3]|5]|3][7|) are tested 
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by measurements of the interaction energy (UJ|5]|P]L70). Anisotropic expansion of strongly- 
interacting Fermi gases Q is analogous to the "elliptic flow" of a quark-gluon plasma ©. 
High temperature superfluidity has been predicted (177117 2117 Jll7 4117 5ll7oT) in strongly-interacting 
Fermi gases, which can be used to test theories of high temperature superconductivity dT7b . 
Microscopic evidence for superfluidity has been obtained by observing the pairing of fermionic 
atoms d7#ll79ll211 . Macroscopic evidence arises in anisotropic expansion (7) and in collective 
excitations (I2/II22II2H) . 

In superconductivity and superfluidity, measurements of the heat capacity have played an 
exceptionally important role in determining phase transitions (24) and in characterizing the na- 
ture of bosonic and fermionic excitations. We report on the measurement of the heat capacity for 
a strongly-interacting Fermi gas of 6 Li atoms, confined in an optical trap. Our experiments (1231) 
examine the fundamental thermodynamics of the gas. 

Thermodynamical properties of the BCS-BEC crossover system are computed (26) using a 
consistent many-body theory (EZIEU based on the conventional mean field state (I2P1) . BCS- 
BEC crossover refers to the smooth change from the Bardeen-Cooper-Schrieffer superfluidity 
of fermions to the Bose-Einstein condensation of dimers, by varying the strength of the pair- 
ing interaction (for example, by tuning a magnetic field). The formalism of Ref. (I76II77II21?!) 
was applied recently (30) to explain radio frequency measurements of the gap (I27J) . The the- 
ory contains two contributions to the entropy and energy arising from fermionic and bosonic 
excitations. The latter are associated principally with excited pairs of fermions (Cooper pairs at 
finite momentum). In this model, there is no direct boson-boson coupling, and fermion-boson 
interactions are responsible for the vanishing of the pair chemical potential n pair in the super- 
fluid regions. The vanishing of n pair implies that, within a trap, the associated low temperature 
power laws in the entropy and energy are the same as those of the homogeneous system d37T) . 
This is to be contrasted with models which involve noninteracting bosons and fermions 021 . 
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Clearly, our BCS-like ground state ansatz will be inapplicable at some point when the fermionic 
degrees of freedom have completely disappeared, and the gas is deep in the BEC regime, where 
the power laws associated with true, interacting bosons are expected (0~7T) . In that case, direct 
inter-boson interactions must be accounted for and they will alter the collective mode behav- 
ior &3l . However, on the basis of collective mode experiments (|2ill22ll23b and their theoretical 
interpretation (1341 1351) . one can argue that the BCS-like ground state appears appropriate in the 
near resonance, unitary regime. The thermodynamic quantities within the trap are computed 
using previously calculated profiles (OoT) of the various energy gaps and the particle density as 
a function of the radius. 

Unlike the weak coupling BCS limit, the pairing gap in the unitary regime is very large. 
Well below the superfluid transition temperature T c , fermions are paired over much of the trap, 
and unpaired fermions are present only at the edges of the trap. These unpaired fermions tend to 
dominate the thermodynamics associated with the fermionic degrees of freedom, and lead to a 
higher (than linear) power law in the temperature (T) dependence of entropy. The contribution 
from finite momentum Cooper pairs leads to a T 3 / 2 dependence of the entropy on temperature. 
Both bosonic and fermionic contributions are important at low T. 

An important feature of these fermionic superfluids is that pair formation occurs at a higher 
temperature T* than the temperature T c where pairs condense. At temperatures T > T*, the 
entropy approaches that of the noninteracting gas. For T c < T < T*, the attraction is strong 
enough to form quasi-bound (or preformed) pairs which are reflected in the thermodynamics. 
At these temperatures, a finite energy, i.e., the pseudogap, is needed to create single fermion 
excitations (I2^II77II7oT) . Interestingly, in the unitary regime, both T* and T c are large frac- 
tions of the Fermi temperature Tp, signifying high temperature pair formation and very high 
temperature superfluidity. 

We prepare a degenerate, unitary Fermi gas comprising a 50-50 mixture of the two lowest 
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spin states of Li atoms near a Feshbach resonance. To cool the gas, we use forced evaporation 
at a bias magnetic field of 840 G in an ultrastable CO2 laser trap (ITll2ll2oT). After cooling well 
into the degenerate regime, energy is precisely added to the trapped gas at fixed atom number, 
as described below. The gas is then allowed to thermalize for 0.1s before being released from 
the trap and imaged at 840 G after 1 ms of expansion to determine the number of atoms and the 
temperature parameter T. For our trap the total number of atoms is iV = 2.2(0.3) x 10 5 . The 
corresponding noninteracting gas Fermi temperature is Tp = {?>N) l ^f\ijj /ks — 2.5 /iK, small 
compared to the final trap depth of Uo/ks = 35 /iK. 

Energy is precisely added to the trapped gas at fixed atom number by releasing the cloud 
from the trap and permitting it to expand for a short time < theat < 460 /is after which the gas 
is recaptured. Even for the strongly-interacting gas, the energy input is well-defined for very low 
initial temperatures, where both the equation of state and the expansion dynamics are known. 
During the times theat used in the experiments, the axial size of the gas changes negligibly, 
while transverse dimensions expand by a factor b±(t heat ). Hence, the mean harmonic trapping 
potential energy (Uho) in eacn of the two transverse directions increases by a factor b\(theat)- 

The initial potential energy is readily determined at zero temperature from the equation of 
state of the gas, (1 + f3)e F (x.) + U HO {x) = ^0 (00), where ep(x) is the local Fermi energy, (3 
is the unitary gas parameter (000|6]0), and /x is the global chemical potential. This equation 
of state is supported by low temperature studies of the breathing mode d2iH23l[J?ll?5t and the 
spatial profiles d7]l6l[5oT). It is equivalent to that of a harmonically trapped noninteracting gas 
of particles with an effective mass (0), which in our notation is m* — m/(l + 0), where m is 
the bare fermion mass. The mean potential energy is half of the total energy, because the gas 
behaves as a harmonic oscillator. As (5 < (00), m* > m, so that the effective oscillation 
frequencies and the chemical potential are simply scaled down, i.e., /xq = ksTp^X + (5 (00). 
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The total energy at zero temperature, which determines the energy scale, is therefore 



E 



4^° = 4 




(1) 



For each direction, the initial potential energy at zero temperature is E /6. Then, the total 
energy of the gas after heating is given by, 



neglecting trap anharmonicity d2oT ). Here, r\ is a correction factor arising from the finite temper- 
ature of the gas prior to the energy input. For the strongly-interacting gas, the initial reduced 
temperature is very low. We assume that it is ~ T = 0.04, where T is measured and calibrated 
as described below. Assuming a Sommerfeld correction then yields r) int ~ l + 27r 2 T 2 /3 ~ 1.01, 
which hardly affects the energy scale. 

A zero temperature strongly-interacting gas expands by a hydrodynamic scale factor b^^theat), 
when released from a harmonic trap (1711771). Heating arises after recapture and subsequent equi- 
libration, but not during expansion. This follows from the lowest T = 0.04, obtained by imaging 
the gas 1 ms after release from the trap. Hence, the temperature change during th eat < 460 /is 
< 1 ms must be very small. 

Thermometry of strongly-interacting Fermi gases is not well understood. By contrast, ther- 
mometry of noninteracting Fermi gases can be simply accomplished by fitting the spatial dis- 
tribution of the cloud (after release and ballistic expansion) with a Thomas-Fermi (T-F) profile, 
which is a function of two parameters. We choose them to be the Fermi radius o x and the 
reduced temperature T/T F . However, this method is only precise at temperatures well be- 
low 0.5 Tp, where a x and T /Tp are determined independently. At higher temperatures, where 
the Maxwell-Boltzmann limit is approached, such a fit determines only the product a 2 . T /Tp. 
We circumvent this problem by determining a x from a low temperature fit, and then hold it 



E{t heat ) =r]E 



+ o b ± (t heat ) , 



(2) 



5 



constant in the fits at all higher temperatures, enabling a one-parameter determination of the 
reduced temperature. 

Spatial profiles of strongly-interacting Fermi gases closely resemble T-F distributions, as ob- 
served experimentally (17117?! and as predicted (136*1) . The profiles of the trapped and released gas 
are related by hydrodynamic scaling to a good approximation. Over a wide temperature range, 
this scaling is consistent with the observed cloud size to ± 2% and is further supported by 
measurements of the breathing frequency, which are within ± 1% of the unitary hydrodynamic 
value (ETJ). Analogous to the noninteracting case, we define an experimental dimensionless 
temperature parameter T, which is determined by fitting the cloud profiles with a T-F distribu- 
tion (Offl) . holding constant the Fermi radius of the interacting gas, a' x . We find experimentally 
that T increases monotonically from the highly degenerate regime to the Maxwell-Boltzmann 
limit. This fitting procedure also leads us to define a natural reduced temperature scale in terms 
of the zero temperature parameters /? and Tp, 

Eq. [21 is consistent with our choice of fixed Fermi radius a' x , i.e., = /i . At high 

temperatures, we must interpret T = T nat , to obtain the correct Maxwell-Boltzmann limit. At 
low temperatures, T ~ T nat yields an estimate of T/Tp which can be further calibrated to the 
theoretical reduced temperature T/Tp by performing the experimental fitting procedure on the 
theoretically generated density profiles (12 71 1261) . 

Preliminary data processing yields normalized, one-dimensional spatial profiles of the atomic 
cloud (26). To determine T over the full temperature range of interest, we employ a fixed expan- 
sion time of 1 ms. We first measure a' x from our lowest temperature data. Then, T is determined 
using the one parameter T-F fit method. This yields T = 0.04 — 2.15 for the strongly-interacting 
gas. 
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The experimental energy scale Eq. [T] and the natural temperature scale Eq. |3] are deter- 
mined by measuring the value of (3. This is accomplished by comparing the measured radius 
of the strongly-interacting gas a' x to the radius for a noninteracting gas d2oT) . We find that 
[3 = —0.49(0.04) (statistical error only) in reasonable agreement with the best current predic- 
tions, where f3 = -0.56 ©, and (3 = -0.545 ©. 

We now apply our energy input and thermometry methods to measure the heat capacity of 
our optically trapped Fermi gas, i.e., for different values of theat, we measure the temperature 
parameter T and calculate the total energy E(theat)/E from Eq.El The time t heat determines 
the energy accurately, as the trap intensity switches in less than 1 /is. We believe that shot-to- 
shot fluctuations in the energy are negligible, based on the small fractional fluctuations in T at 
low temperatures, where the heat capacity is expected to be very small. To obtain high resolution 
data, 30-40 different heating times theat are chosen. The data for each of these heating times are 
acquired in a random order to minimize systematic error. Ten complete runs are taken through 
the entire random sequence. 

We first measure the heat capacity for a noninteracting Fermi gas (|2ill26l) . where the scat- 
tering length a is zero. This occurs near 526 G. Fig. 1 shows the data (green dots) which 
represent the calculated E(t heat )/E versus the measured value of T, for each t heat . For com- 
parison, predictions for a noninteracting, trapped Fermi gas, Eid ea i{T) / Eid ea i(0) are shown 
as the black curve, where T = T/T F in this case. Here, the chemical potential and energy 
are calculated using a finite temperature Fermi distribution and the density of states for the 
trapped gas. Throughout, we use the density of states for a realistic Gaussian potential well, 
JJ{r) = U [l - e - m ^ 2r2 / 2C/ °] with U = lA.Qk B T F , rather than the harmonic oscillator ap- 
proximation. This model is in very good agreement with the noninteracting gas data at all 
temperatures. 

For the strongly-interacting gas at 840 G, Fig. 1 (blue diamonds), the gas is cooled to T = 
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0.04 and then heated. Note that the temperature parameter T varies by a factor of 50 and the 
total energy by a factor of 10. For comparison, we show the theoretical results for the unitary 
case as the red curve. Here the horizontal axis for the theory is obtained using the approximation 
T ~ T nat via Eq.[3j On a large scale plot, the data for the strongly-interacting and noninteracting 
gases appear quite similar, although there are important differences at low temperature. 

A striking result is observed by plotting the low temperature data of Fig. 1 on an expanded 
scale (I25II26I) . This reveals a transition in the heat capacity which is made evident by plotting the 
data for the strongly-interacting gas on a log — log scale as in Fig. 2. The transition is apparent 
in the raw temperature data d2"5ll2oT) . and is strongly suggestive of the onset of superfluidity. 
Note that the observed spatial profiles of the gas vary smoothly and are closely approximated 
by T-F shapes in the transition region. Fig. 2 shows the transition after converting the empirical 
temperature T to theoretical T/T F units. 

The empirical temperature is calibrated to enable precise comparison between the theory 
and the experimental data. For the calibration, we subject the theoretically derived density pro- 
files 13611271) to the same one-dimensional T-F fitting procedure as used in the experiments. One 
dimensional density distributions are obtained by integrating over two of the three dimensions 
of the predicted spatial profiles, which are determined for a spherically symmetric trap. Our 
results for this temperature calibration are shown in the inset to Fig. 2. This calibration pro- 
vides a mapping between the experimental reduced temperature a/1 + P T and the theoretical 
temperature T/T F . We find that T = T nat is a very good approximation above T c . Such scaling 
may be a manifestation of universal thermodynamics (4). The difference between T and T nat is 
significant only below the superfluid transition T c and is therefore negligible in the large scale 
plot of Fig. 1 over a broad temperature range. However, below T c the fits to the theoretical 
profiles yield a value of a/1 + (3 T which is lower than the theoretical value of T/Tp. This is a 
consequence of condensate effects (26). 



8 



Fig. 2 shows that above a certain temperature T c , the strongly-interacting data nearly overlap 
that of the noninteracting gas, and exhibit a power law fit E/Eq — 1 = 4.98(T/Tp) 143 . Below 
T c , the data deviate significantly from noninteracting Fermi gas behavior, and are well fit by 
E/Eq — 1 = 97.3(T/Tp) 3 73 (dashed curve). From the intersection point of these power law 
fits, we estimate T c /Tp = 0.27(.02) (statistical error only). This is very close to our theoretical 
value TJT F = 0.29. 

The fractional change in the heat capacity C is estimated from the slope change in the fits 
to the calibrated data. In that case, the relative specific heat jump (C< — C > )/C > ~ 1.51(0.05) 
(statistical error only), where > (<) denotes above (below) T c . This is close to the value (1.43) 
for an s-wave BCS superconductor in a homogeneous case, although one expects pre-formed 
pairs, i.e., pseudogap effects, to modify the discontinuity somewhat 

In Fig. 2 and Fig. 3, the theory is compared to the calibrated data after very slightly de- 
tuning the magnetic field in the model away from resonance, so that the predicted unitary gas 
parameter (3 has the same value as measured. This small detuning, [kpa)^ 1 = 0.11, where 
kp = \j2mkBTp/h 2 , is reasonable given the broad Feshbach resonance (OPI) in 6 Li. 

Finally, Fig. 3 presents an expanded view of the low temperature region. Here, the ex- 
perimental unitary data is calibrated and replotted in the more conventional theoretical units, 
Ep = ksTp and Tp. The agreement between theory and experiment is very good. In the pres- 
ence of a pseudogap, a more elaborate treatment (28) of the pseudogap self-energy, which takes 
into account spectral broadening, will be needed in order to calculate accurately the specific 
heat jump. 

If one extends the temperature range in Fig. 3 to high T we find that both the unitary and 
noninteracting cases coincide above a characteristic temperature, T*, although below T c they 
start out with different power laws (as shown in Fig. 2). In general, we find that agreement 
between theory and experiment is very good over the full temperature range for which the data 
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are taken. The observation that the interacting and noninteracting curves do not precisely co- 
incide until temperatures significantly above T c is consistent with (although it does not prove) 
the existence of a pseudogap and with onset temperature from the figure T* « 2T C . Related 
signatures of pseudogap effects are also seen in the thermodynamics of high temperature super- 
conductors (071) . 
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Figure 1 : Total energy versus temperature. For each heating time theat, the temperature parame- 
ter T is measured from the cloud profile, and the total energy E(theat) is calculated from Eq. © 
in units of the ground state energy E . Green circles: noninteracting Fermi gas data; Blue di- 
amonds: strongly-interacting Fermi gas data. Black curve: predicted energy versus reduced 
temperature for a noninteracting, trapped Fermi gas, E idea i(T) / E idea i(0); Red curve: predicted 
energy versus T for the unitary case. No temperature calibration is applied since T « T nat 
over the broad temperature range shown. Note that the lowest temperature point (blue square) 
is constrained to lie on the black curve. 
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Figure 2: Energy input versus temperature from Fig. 1 after temperature calibration on a log — 
log scale. The strongly-interacting Fermi gas shows a transition in behavior near T/T F = 0.27. 
Green circles: noninteracting Fermi gas data; Blue diamonds: strongly-interacting Fermi gas 
data; Red (Black) curve: prediction for a unitary (noninteracting), Fermi gas in a Gaussian trap 
as in experiment; Black dashed line: best fit power law 97.3 (T/Tp) 3 73 to the unitary data for 
T/T F < 0.27. The inset shows the calibration curve, which has been applied to the unitary data 
(blue diamonds). The red dashed line in the inset represents the diagonal, T/Tp = y/T+~j3T. 
Hereto = E(T = 0). 
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Figure 3: Low temperature comparison of present theory (red, black curves) and experiments 
(symbols) in terms of E/Ep (Ep = ksTp) per atom as a function of T/Tp, for both unitary 
and noninteracting gases in a Gaussian trap. The fact that the two experimental (and the two 
theoretical) curves do not merge until higher T* > T c is consistent with the presence of a 
pseudogap. 
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Supporting Online Material 
Computation of Thermodynamical Quantities 

The theoretical community is in the midst of unraveling the nature of resonantly interacting 
fermionic superfluids (,5|7J $2\ $3\ &\ $5\ $6\ 43 $£\ -® with particular emphasis on the strongly 
interacting Fermi gas (Silk. In the BCS-BEC crossover picture (S ll), the strongly interacting 
Fermi gas is intermediate between the weak coupling BCS and BEC limits. In addressing the 
nature of the excitations from the conventional mean field or BCS-like ground state (S12), our 
theoretical calculations help to provide a theoretical calibration of the experimental thermome- 
try, and elucidate the thermodynamics. 

Without doing any calculations one can anticipate a number of features of thermodynam- 
ics in the crossover scenario. The excitations are entirely bosonic in the BEC regime, exclu- 
sively fermionic in the BCS regime, and in between both types of excitation are present. In 
the so-called one-channel problem the "bosons" correspond to noncondensed Cooper pairs, 
whereas in two-channel models, these Cooper pairs are strongly hybridized with the molecular 
bosons of the closed channel, singlet state. Below T c the presence of the condensate leads to 
a single-branch bosonic excitation spectrum which, at intermediate coupling, is predominantly 
composed of large Cooper pairs. These latter bosons lead to a pseudogap ( £771 above T c . 
Within the conventional mean field ground state, and over the entire crossover regime (S14) 
below T c , the bosons with effective mass M* have dispersion tt q = h 2 q 2 /2M*. This form for 
the dispersion reflects the absence of direct boson-boson interactions. In the extreme BEC limit, 
when the fermionic degrees of freedom become irrelevant, direct inter-boson interactions must 
be accounted for. While our focus in this paper is on the unitary case, when we refer to "BEC" 
we restrict our attention to the near-unitary BEC regime. 

As long as the attractive interactions are stronger than those of the BCS regime, these non- 
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condensed pairs must show up in thermodynamics, as must the pseudogap in the fermionic 
spectrum. These are two sides of the same coin. Below T c , the fermionic excitations have dis- 
persion Ek = a/ (ek — /i) 2 + A 2 , where ek = h 2 k 2 /2m and jj, are the atomic kinetic energy 
and fermionic chemical potential, respectively. That this excitation gap A is non-zero at T c 
in the Bogoliubov quasi-particle spectrum E^, differentiates the present approach (S14) from 
all other schemes which address BCS-BEC crossover at finite T. The bosons, by contrast, are 
gapless in the superfluid phase, due to their vanishing chemical potential. Within a trap, and 
in the fermionic regime (for which ji > 0), the fermionic component will have a strong spatial 
inhomogeneity via the spatial variation of the gap. Thus, in contrast to the homogeneous case, 
fermions on the edge of the trap, which have relatively small or vanishing excitation gaps A, 
will contribute power law dependences to the thermodynamics. 

Starting at a magnetic field well above a Feshbach resonance, by decreasing the magnetic 
field, we tune from the BCS-like regime towards unitarity at resonance. We first consider low T 
where fermions become paired over much of the trap. The unpaired fermions at the edge tend 
to dominate the thermodynamics associated with the fermionic degrees of freedom, and lead 
to a higher (than linear) power law in the T dependence of the entropy. The contribution from 
excited pairs of fermions is associated with a T 3 / 2 dependence of entropy on temperature which 
dominates for temperatures T/Tp < 0.05 or T/T c < 0.2. In general, the overall exponent 
of the low T power law varies with magnetic field, depending on the magnitude of the gap 
and temperature, as well as the relative weight of fermionic and bosonic contributions. In the 
superfluid phase, at all but the lowest temperatures, the fermions and bosons combine to yield 
S oc T 2 precisely at resonance ((fcpa) -1 = 0). For the near-unitary case investigated in the 
paper {(kpa)' 1 = 0.11), we have S oc T 19 . 

Because our calculations C3751) are based on the standard mean field ground state G572J, we 
differ from other work ($Z\ < 3i6l) at finite temperatures. Elsewhere ST3\ STh we have char- 
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acterized in quantitative detail the characteristic gap A and pseudogap A pg energy scales. The 
pseudogap (which is to be associated with a hybridized mix of noncondensed fermion pairs and 
molecular bosons) and the superfluid condensate (sc) called A sc , add in quadrature to determine 
the fermionic excitation spectrum: A 2 (T) = A 2 C (T) + A 2 pg (T). Our past work $B\ 477ti 
has primarily focussed below T c . Here we extend these results, albeit approximately, above T c . 
Our formalism has been applied below T c with some success in Ref. (.5© to measurements of 
the pairing gap in RF spectroscopy. A more precise, but numerically more complex method for 
addressing the normal state was given in Ref. (S18). 

After including the trap potential U(r) and internal binding energy of the bosons, the local 
energy density can be decomposed into fermionic (Ef) and bosonic (Ef,) contributions and 
directly computed as follows 

E = /m(r) +E f + E b , 
E f = J2( iaj n + e k -v(r))G(K) 

K 

= ^[2£ k /(£ k ) - (E k - e k + /i(r))] + A 2 X (0) , 

k 

Eb = ^^(fiq — Hboson) H^q — ^boson) , (SI) 

q 

where fi(r) = fi—U(r), n(r) is the local density, tu n = (2n+ l)n ksT is the fermionic Matsubara 
frequency, G(K) is the renormalized fermionic Green's function with four-momentum K = 
(iu n , k), b(x) and f(x) are the Bose and Fermi distribution functions, respectively. The pair 
susceptibility x(0), at zero frequency and zero momentum, is given by 

k k 

and the bosonic chemical potential /i& oson is zero in the superfluid phase. 

Unlike the situation in condensed matter systems, for these ultracold gases, thermometry 
is less straightforward. Experimentally, temperature is determined from the spatial profiles 
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of the cold gas, either in the trap, or following expansion. For weakly interacting Bose and 
Fermi gases, where the theoretical density is well understood, this procedure is straightforward. 
However, for a strongly interacting gas, the spatial profile has not been understood until re- 
cently (STh. For this reason, the temperature is often measured on either side far away from 
the Feshbach resonance, where the scattering length is small. A strongly interacting sample in 
the unitary regime is then prepared by an adiabatic change of the magnetic field. 

More specifically, in the BCS or weak attraction regime, temperature is determined by fitting 
the spatial (or momentum distribution) profiles to those of a non-interacting Fermi gas In 
the opposite BEC regime, temperature can be deduced by fitting the Gaussian wings of density 
profiles or determining condensate fractions ( $2U\ $2l\ . Thus, it is convenient to describe a 
given intermediate regime which is accessed adiabatically, by giving the initial temperature 
at either endpoint. In order to determine this adiabatically accessed temperature, one needs 
precise knowledge of the entropy S as a function of T and magnetic field from BCS to BEC. 
The entropy S can be calculated directly CSl75t as a sum of fermionic and bosonic contributions 
based on the two types of excitations. Equivalently, one can also calculate the entropy from the 
energy, S = £ f§. 

In the strongly interacting regime, one can measure an empirical temperature T by fitting a 
T-F density profile directly to the spatial distribution, as done in this paper. In the following, we 
describe a temperature calibration method which relates the measured empirical temperature T 
to the theoretical value of T/T F . 

Calibration of Experimental Temperature Scale 

In order to obtain a temperature calibration curve for the experiments (inset, Fig. 2 main text) we 
note that our theoretically generated profiles yield very good agreement with the Thomas-Fermi 
functional form (S17) for the normal and superfluid states. However, there are slight systematic 
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deviations from this form in the superfluid phase. Below T c the profiles contain the superfluid 
condensate as well as non-condensed pairs along with excited fermions. Although our profiles 
are generated for an isotropic trap, it can easily be shown that trap anisotropy is not relevant for 
thermodynamic quantities. Because they involve integrals over the entire trap, the calculations 
can be mapped onto an equivalent isotropic system. 

Our theoretical profiles are generated for given reduced temperatures T/T F . If one applies 
the experimental procedure to these theoretical profiles one can deduce the parameter a/1 + ]3 T 
for each T/Tp. Theoretically, then, it is possible to relate these two temperature scales. This is 
summarized by the calibration curve in the inset to Figure 2. 

Quite remarkably, it can be seen from this inset that the experimental T-F fitting procedure 
yields the precise theoretical temperature in the normal state. This applies even below the pseu- 
dogap onset temperature T*, since the non-condensed pairs and the fermions both are thermally 
distributed. However, in the superfluid phase, the parameter a/1 + (5 T systematically under- 
estimates the temperature, because of the presence of a condensate. One can understand this 
effect as arising principally from the fact that the region of the trap occupied by the conden- 
sate is at the center and decreases in radius as temperature is increased, until it vanishes at T c . 
This prevents the profile from expanding with temperature as rapidly as for the non-interacting 
fermions of strict T-F theory. Hence, one infers an apparently lower temperature. As T/Tp 
approaches zero, the parameter a/1 + ]3 T must approach zero as well. 

Experimental Methods and Empirical Thermometry 

Preparation of the strongly interacting Fermi gas is described in the main text and the details 
can be found elsewhere ^21 $23. 

Preparation of degenerate, noninteracting Fermi gases follows a similar series of steps. As 
described previously Cfl2"21). 23 s of forced evaporation at 300 G brings the temperature of the 
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gas to T = 0.24, the lowest temperature we can achieve in this case. The gas is then heated as 
described in the main text. Finally, the gas is released and imaged at 526 G to determine the 
number of atoms and the temperature. Temperatures T between 0.24 and 1.23 are obtained for 
the noninteracting gas. 

All heating and release for time of flight measurements are conducted at 4.6% of the full 
trap depth. At this depth, the measured trap frequencies, corrected for anharmonicity, are u± = 
y/UJJJ^ = 2nx 1696(10) Hz and uj z = 2n x 72(5) Hz, so that u> = (u} x u}yU z ) 1/3 = 2nx 592(14) 
Hz is the mean oscillation frequency. 

For both the interacting and noninteracting samples, the column density is obtained by 
absorption imaging of the expanded cloud after 1 ms time of flight, using a two-level state- 
selective cycling transition ( Sffift SU2~[ . In the measurements, we take optical saturation into 
account exactly and arrange to have very small optical pumping out of the two-level system. 
The resulting absorption image of the cloud can then be analyzed to determine the temperature 
of the sample. 

Anharmonic Corrections to the Energy Input 

Eq. 2 of the main text does not include corrections to the energy input which arise from anhar- 
monicity in the gaussian beam trapping potential. In general, after the cloud expands for a time 
theat, the energy changes when the trapping potential U (x) is abruptly restored, 



Here n(x, theat) 0^o( x )) is the density of the expanded (trapped) cloud, where no(x) is a zero 
temperature T-F profile, as noted in the main text. A scale transformation (SI (ft SHty relates 
n (x, theat ) to no (x) . Using this result, we obtain Eq. 2 of the main text as well as the anharmonic 
correction AE arising for a gaussian beam trapping potential. For a cylindrically symmetric 




(S3) 
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trap, we obtain, 



AE 



[2bi(t) + bl(t)-3] + 



[ 46i (t) + 2b\ (t) + 361 0) - 9 ] ■ (S4) 



30 a 



o 



360 Ul 



Note that for our experiments, we assume a gaussian beam potential with three different di- 
mensions. These corrections are most significant for the largest values of theat, since the largest 
contribution to the energy change arises from atoms at the edges of the cloud. 

Energy Input for Noninteracting Samples 

Although the interacting and noninteracting samples are heated in the same fashion, there are 
a few differences in the way the energy input is calculated. In the noninteracting case, the 
correction factor in Eq. 2 of the main text, r] nonint , is determined at the lowest temperature T = 
0.24 from the energy for an ideal Fermi gas. Furthermore, whereas the strongly interacting gas 
expands hydrodynamically, expansion of the noninteracting gas is ballistic so that b±(th ea t) = 

bf (theat) = a/ 1 + Ntatf- 

Determination of (3 

We determine (3 by comparing the measured Fermi radius for the strongly interacting sample 
a x to the calculated radius for a noninteracting gas o x confined in the same potential. The 
relation is given by a x — a x (l + /3) 1 / 4 ($2~5\i. where a x = ^j2k B T F / (Muj x ) is the radius for a 
noninteracting gas. We obtain o x = 1.065 (iV/2) 1//6 fxm for our trap parameters. This calculated 
radius is consistent with the value measured for noninteracting samples at 526 G in our trap. To 
determine a' x , we measure the size of the cloud after 1 ms of expansion, and scale it down by the 
known hydrodynamic expansion factor of 6^(1 ms) = 13.3 (. 37(31 We then determine the 
Fermi radius a' x = 11.98 (N/2) 1 / 6 /im/13.3 = 0.901(0.021) (A^/2) 1 / 6 /im. With these results, 
we obtain [3 = —0.49(0.04) (statistical error only). 
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Observed Transition in Energy versus Empirical Temperature T 

For the strongly interacting Fermi gas, without calibrating the empirical temperature scale, we 
observe a transition between two patterns of behavior at T = 0.33 (^2oT): For T = 0.33 — 2.15, 
we find that the energy closely corresponds to that of a trapped Fermi gas of noninteracting 
atoms with the mass scaled by 1/(1 + 0). At temperatures between T = 0.04 — 0.33, the energy 
scales as T 2 53 , significantly deviating from ideal gas behavior as can be seen in Fig. [ST| The 
transition between two power laws is evident in the slope change of the log — log plot of Fig.lS21 
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Figure SI: Strongly-interacting Fermi gas below the transition temperature. E/E versus un- 
calibrated empirical temperature T on a linear scale. Orange line, best fit power law 9.8 T 2 53 . 
Black curve: Predicted E/E for an ideal Fermi gas as a function of T = T/Tp. Note the 
lowest temperature point (blue square) is not included in the fits: It is constrained to lie on the 
black curve by our choice of rj int = 1.01 in Eq. 2 of the main text. 
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Figure S2: Energy input versus uncalibrated temperature T on a log — log scale. The strongly in- 
teracting Fermi gas shows a transition in behavior near T = 0.33. Green circles: noninteracting 
Fermi gas data; Blue diamonds: strongly interacting Fermi gas data. Black curve, prediction for 
a noninteracting, trapped Fermi gas. Orange line, best fit power law 9.8 T 2 53 . Note the lowest 
temperature point (blue square) is not included in the fits, as it is constrained to lie on the black 
curve. 
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